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Abstract 

A model is presented in which supernova-driven winds from early 
galaxies reionize the intergalactic medium by z = 5. This scenario can 
explain the observed absence of a Gunn-Peterson trough in the spectra 
of high-redshift quasars providing that the bulk of these early galaxies 
are quite small, no more massive than about 10 8 M Q . It also predicts 
that most of the IGM was enriched to at least 10% of current metal 
content by z — 5 and perhaps as early as z = 15. The existence of such 
early mini-galaxies violates no spectral constraints and is consistent 
with a pure CDM model with b < 2. Since the final radius of a typical 
ionized bubble is only around 100 kpc, the induced modification of the 
galaxy autocorrelation function is negligible, as is the induced angular 
smoothing of the CBR. Some of the gas swept up by shells may be 
observable as pressure-supported Lyman-alpha forest clouds. 
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1 Introduction 



The absence of a Gunn-Peterson trough in the spectra of high-redshift 
quasars has provided strong evidence for the intergalactic medium (IGM) 
being highly ionized as early as z = 4 (Gunn & Peterson 1965; Steidel and 
Sargent 1987; Webb et al. 1992). The hypothesis that photoionization of the 
IGM by quasars could account for this ionization (Arons & McCray 1969; 
Bergeron &: Salpeter 1970; Sherman 1980) has been challenged (Shapiro 
1986; Shapiro & Giroux 1987; Miralda-Escude k Ostriker 1990). Other 
studies have maintained that photoionization by quasars (Donahue & Shull 
1987) or active galactic nuclei (Teresawa 1992) may nonetheless be sufficient. 
However, in view of the large uncertainties in crucial parameters such as ion- 
izing fluxes, the issue of what reionized the IGM must still be considered 
open. 

In comparing the Gunn-Peterson constraints with our work in papers 3 
and 4, the crucial difference is the degree of ionization required. To affect 
the CBR, it does not really matter whether the ionization fraction x is 90% 
or 99.999%, as this makes only a 10% difference in the optical depth r. The 
Gunn-Peterson limits constrain not x but (1 — x), the neutral fraction. Thus 
in this context, the difference between 90% and 99.999% ionization is four 
orders of magnitude. In papers 3 and 4, we found that photoionization by 
early galaxies could easily reionize the IGM by a redshift z = 5, but the 
issue here is whether photoionization alone can provide the extremely high 
ionization fraction required to pass the Gunn- Peterson test. 

In this paper, we investigate an alternative reionization scenario, which 
produces considerably higher IGM temperatures than those attained by the 
photoionization models in previous papers. Supernova driven winds from 
luminous galaxies have long been conjectured to be an important ionization 
source for the IGM (Schwartz et al. 1975; Ikeuchi & Ostriker 1986; Carlberg 
& Couchman 1989). Cold dark matter (CDM)-based models of structure 
formation (Blumenthal et al. 1984; Efstathiou et al. 1985) predict the for- 
mation of gravitationally bound objects of mass as small as 10 7 M & in large 
numbers before z = 5. Recent work (Blanchard et al. 1992) indicates that 
such objects can cool rapidly and presumably fragment into stars. These 
early mini-galaxies would be expected to release great amounts of kinetic 
energy into the surrounding IGM, thereby creating large, fairly spherical 
voids filled with thin, hot, ionized plasma. We analyze the effect of expand- 
ing bubbles driven by supernova winds from early mini-galaxies, and show 
that this mechanism of distributing energy can indeed provide the required 
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ionization without violating any of the current spectral constraints. 
In Section ^, we will treat the expansion of a shell in a uniform, cold and 
neutral IGM. As these bubbles become larger and more numerous and fill 
most of space, this obviously becomes a very poor model of the IGM. In 
Section ^ we estimate bulk properties of this new processed IGM such as 
temperature, density and ionization. 

2 The Explosion Model 

Since the pioneering work on spherically symmetric explosions by Sedov 
(1959), a profusion of analytic solutions have been given by numerous au- 
thors for models of ever-increasing complexity (Cox & Smith 1974; McKee 
& Ostriker 1977; Weaver et al. 1977; McCray & Snow 1979; Bruhweiler et 
al. 1980; Tomisaka et al. 1980; McCray & Kafatos 1987; Ostriker & McKee 
1988). Most of these models pertain to bubbles in the interstellar medium 
of a galaxy, where the expansion of the universe can be ignored. Ostriker 
& McKee have given asymptotic self-similarity solutions that incorporate 
this latter complication, but unfortunately they are not sufficiently accurate 
for our needs. The reason is that since neither energy nor momentum is 
conserved in the regime before the shell becomes self-similar, there is no 
accurate way to normalize the self-similar solution using the initial data. 

Let pt, and pd denote the average densities of baryonic and non-baryonic 
matter in the universe. We will assume that all baryons are in diffuse for 
early on, so that pb is also the density of the IGM. We will write pb = QbPc 
and pd = Qdp c , where the critical density p c = 3H 2 /8ttG. We will use a 
three-phase model for the expanding bubbles: 

• a dense, fairly cool spherical shell of outer radius R and thickness R5, 
containing a fraction (1 — f m ) of the total baryonic mass enclosed, 

• uniform neutral ambient intergalactic medium (IGM) of density Pb+pd 
and zero pressure outside, 

• a hot, thin, isothermal plasma of pressure p and temperature T inside 
the shell. 

The shell is driven outwards by the pressure of the hot interior plasma 
but slowed by the IGM and by gravity. The plasma is heated by kinetic 
energy from supernova explosions and collisions with IGM and cooled by 
bremsstrahlung and Compton drag against the cosmic background radiation. 
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2.1 The expanding shell 

We assume that the expanding shell sweeps up almost all baryonic IGM that 
it encounters and loses only a small fraction of it through evaporation into 
the interior, so that its total mass is given by m(t) = |vri?(t) 3 (l — f m )Pb, 
where the constant f m <C 1. Since Pb/Pb = —3H for any cosmological model, 
we get 

— = (R 3 pb) 1 4" ( R3 Pb) = 3 ( ^ - H ) if ^ > H, zero otherwise. 
m V / at V / \it / it 

(The shell will acquire new mass when it is expanding faster than the Hubble 
flow, and will never lose mass.) It turns out that the Hubble flow catches up 
with the shell only as t — ► oo, so we will always have R > HR and rh > 0. 

When new mass is swept up, it must be accelerated from the velocity HR 
to R, so the shell experiences a net braking force (R — HR)m. The interior 
pressure p drives the shell outward with a force pA = A7rR 2 p = 3mp/pbR 
in the thin shell approximation 5, f m <C 1. Finally there is a gravitational 
braking force, which in the thin-shell approximation (Ostriker & McKee 
1988) gives the deceleration ^irGR(pd + \pb)- Adding these three force 
terms, the radial equation of motion becomes 

n b H 2 R R\ ) V 2 V 2 K ' 

2.2 The interior plasma 

The equation of state for the plasma in the hot interior gives the thermal 
energy 

E t = \pV = 2irpR\ (2) 

and energy conservation for the interior yields 

E t = L- pdV/dt = L — 4irpR 2 R, (3) 

where the luminosity L incorporates all sources of heating and cooling of 
the plasma. We will consider five contributions to L and write 

L = Lsn ~ L com p — Lbrerns -t'ion ~\~ Ldi ss , 

where L sn is the energy injection from supernova explosions, L comp the cool- 
ing by Compton drag against the CBR, Lb rem s the cooling by bremsstrahlung, 
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Li on the cooling by ionization of neutral IGM and Ldiss the heating from 
collisions between the shell and the IGM. 

In stellar burning from zero to solar metallicity, the mass fraction 0.02 x 
0.007 is released, mostly as radiation. Due to low cross-sections, only a 
negligible fraction of this radiation will contribute towards heating the gas, 
so we will only be interested in the energy that is released in kinetic form. 
From empirical observations of active galactic winds (Heckman 1990) about 
2% of the total luminosity from a galaxy is mechanical. Another empirical 
observation is that for a solar neighborhood initial stellar mass function, 
one has roughly one supernova for every 15OM of baryons that form stars, 
with a typical kinetic energy output of 10 51 ergs per explosion. Both of 
these observations lead to the same estimate 

L sn = — 0(tbum — t)m 1.2Lo——9(ti )Urn — t), 

tburn M Q 

where the efficiency f sn ~ 4 x 10~ 6 and where we have assumed that the 
energy is released at a constant rate during a period i& urn « 5 x 10 7 years. 

Now let us examine cooling. The interior baryon density is pi = pbf m /(l— 
<5) 3 whereas the shell density is p s = pb(l — f m )/(l — (1 — <5) 3 ) « p&/3<5 if 
f m ,5 <C 1. Compton drag against the microwave background radiation 
causes energy loss at a rate (Kompaneets 1957) 

4tt 2 , , / kT e \ /£;T 7 \ 4 
L comp = — (a t cn e ) *cV, (4) 

where at is the Thomson cross section, V = |7ri? 3 , and T e = T, the temper- 
ature of the interior plasma, which is given by 

Et = (- + -) kT^V. 

(We will assume almost complete ionization and low metallicity, so that 
n e ~ fm n b-) Using equation (g), we see that Compton drag causes cooling 
on a timescale 

Ef 45 / tic \ 4 m P , 9 , . , 

T J - = T2\Vr) -| «2 x 10 12 years x(l + z)~ 4 , 

Lcomp ^ \ K-T 7 J (T t Ti 

that is, it becomes important only at high redshifts. It turns out that 
rems ^ L com p in our regime of interest, so we will simply make the ap- 
proximation L\) Tems ~ 0. Assuming that the ambient IGM is completely 
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neutral, the power required to ionize the hydrogen entering the interior is 
simply 

Lion = fmribEo x 4vri? 2 R - HR , 

where E ~ 13.6 eV. 

The equation of motion ([l]) assumes that the collisions between the ex- 
panding shell and the ambient IGM are perfectly inelastic. The kinetic 
energy dissipated has one of three fates: It may 

(a) radiate away in shock cooling processes, 

(b) ionize the swept up IGM, or 

(c) heat the shell and by conduction heat the interior plasma. 

Let fd denote the fraction that is reinjected into the interior plasma through 
processes (b) and (c). This is one of the major uncertainties of the model. 
Now a straightforward kinematic calculation of the kinetic energy loss per 
unit time gives 

Ldiss = /rf^n ~ HRj ■ 

Making accurate estimates of fd is difficult, so we simply use the two extreme 
cases fd = and fd = 1 in the simulations. Perhaps surprisingly, the results 
will be seen to be relatively independent of the choice of fd- 

2.3 Solutions to the equations 

Combining (|2|) and (|3|) leaves us with 

* = 2^3 " ( 5 ) 

The system of equations (|l|) and @ reduces to that derived by Weaver et 
al. (1977) in the special case where L(t) is constant and the expansion of 
the universe is ignored. 

Let us define dimensionless variables as follows: 

t= t/U, U = | J ff " 1 ( 1 + z *) _3/2 
T] = H/H*, = 

£ = L/L„ U = f sn M b c 2 t^ rn « 1.2 x 1O 5 L x M 5 
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e = E/E+, E* = LJtorn « 7.2 x 10 53 ergs x M f 
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r 



R*= Ll /S G^ B U ps0.13M P cx h- 1 {l + z^-'^Ml 1 ^ 

p/p*, p* = L* /5 G~ 3 / 5 C 2 ~ 1.4 x 1(T 16 Pa x /i 2 (l + z*) 3 M 5 2/5 
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Here we have taken h = 0.5 and denned M 5 M fe /10 5 M . If = Q b + 
tld=l, then is the age of the universe at the redshift z* when the shell 
begins its expansion, i.e. the Big Bang occurred at r = — 1 and the shell 
starts expanding at r = 0. For this simple case, we have r\ = (1 + t) _1 = 

(l + Z,)- 3 / 2 (l + z) 3 / 2 . 

Equations (|l|) and @ now become 



r"ir) = ^V(r)-^ ~ 3 (l - |^) 2 - (§O d + ±0 6 ) ,(r) 2 r(r) 

(6) 

Here £ = £comp ^ion "I - i-dissi where 

£ comp ps 0.017/i- 1 (l + z,)- 3 / 2 (l + z) 4 r 3 g, 

£ ion ps 2.2f m n b M~ 2/5 x r/ 2 r 2 (r' - fryr) , and 

4iss = |/<A> X (V) 2 - lw) • 

In computing £ CO mp, we have taken T 7 o = 2.74K. The interior temperature, 
the thermal energy and the kinetic energy are given by 

T ps 4.5 x W 5 K x 



,5 v . ^5 /5 



f m Q, b 7](t) 2 ' 

E t = r°q, 

7~burn 

mk 2 /2 n b 2 3 /2 

L*t burn 9r b 

urn 

The solution to the system (|6|) evolves through three qualitatively dif- 
ferent regimes: t <C 1, t ps 1 and r> 1. 

a) In the limit of small times t< 1, gravity and Hubble flow are negligible 
and we obtain the asymptotic power law solution 



r(T) = ar 



3/5 



z (t) = 6r" 4 / 5 
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/ 375M A 1/5 , , 7fi 6 2 
Wherea "l77^7^j and 6 =150^' ^ 

as may be verified by direct substitution. This solution reduces to 
that found by Weaver et al. in the special case fd = 0. Since the total 
energy injected is simply £j„ = T/Tb urn , this gives 

e t 35 e k 15 

and — 



E in e m 77- 27/ d 

for small r. Hence even though Et + Ek = £ m only for the most opti- 
mistic case fa = 1, we see that no more than 1 — 35 ,+ 15 « 25% of the 
injected energy is lost as radiation even in the worst case fd = 0. 

b) The behavior in the intermediate regime is a complicated interplay be- 

tween several different effects: 

1. After approximately 5 x 10 7 years, the supernova explosions cease, 
which slows the expansion. In this pressure-driven snowplow 
phase, we would asymptotically have R oc t 2 / 7 , t 4 / 13 , i 1 / 3 , i 4 / 11 
or i 2 / 5 if there were no gravity, no Hubble flow and no cooling 
with fd = 0, |, |, | or 1, respectively. 

2. Cooling (and pdV) work reduces the pressure and the thermal 
energy to virtually zero, which slows the expansion. With zero 
pressure, we would approach the momentum-conserving snow- 
plow solution R oc t 1 / 4 if there were no gravity and no Hubble 
flow. 

3. The density of the IGM drops and the IGM already has an out- 
ward Hubble velocity before it gets swept up, which boosts the 
expansion and adds kinetic energy to the shell. 

4. Gravity slows the expansion. 

5. Dark matter that has been accelerated outward by the shell catches 
up with it again and speeds up the expansion. (This last effect has 
been neglected in the equations above, since it generally happens 
too late to be of importance for our purposes.) 

c) As t — > oo, the shell gets frozen into the Hubble flow, i.e. R oc i 2 / 3 

if Q = 1. An approximate analytic solution for r S> 1 is given by 
Ostriker& McKee (1988), but since neither energy nor momentum is 
conserved in the intermediate regime, there is no simple way to connect 
this solution with the short-time solution above. 
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Numerical solutions for the comoving radius (1 + z)R are plotted in Figure [l] 
for different values of z* and fa- The asymptotic solution (0) has been used 
to generate initial data at r = 0.01 for the numerical integration. In this 
Figure, we have truncated R when the interior temperature drops below 
15,000K, after which newly swept up IGM fails to become ionized. Figures ^, 
H and H] show what becomes of the injected energy for different parameter 
values. Note that the relative fractions are approximately constant early 
on, while the supernovae inject energy, in accordance with the asymptotic 
solution (0). The reason that the total energy exceeds 100% of the input 
is that the shell gobbles up kinetic energy from swept-up IGM that already 
has an outward Hubble velocity. 



3 Cosmological Consequences 

Once the expanding bubbles discussed in the previous section have pene- 
trated most of space, the IGM will presumably have a frothy character on 
scales of a few 100 kpc, containing thick and fairly cool shell fragments 
separated by large, hot, thin and ionized regions that used to be bubble 
interiors. 



In Section 3.1, we calculate at what point the IGM becomes frothy, more 
specifically what fraction of space is covered by expanding shells at each z. 
In |3.2| we discuss the resulting enrichment of the IGM with heavy elements. 
In |3.3| the thermal history of the IGM after this epoch is treated. Finally, 
in |3.4] the residual ionization is computed, given this thermal history, and 
we discuss the circumstances under which the Gunn-Peterson constraint is 
satisfied. 

3.1 IGM porosity 

Assuming the standard PS theory of structure formation (Press & Schechter 
1974), the fraction of all mass that has formed gravitationally bound objects 
of total (baryonic and non-baryonic) mass greater than M at redshift z is 



1 — erf 



V2a(M) 



where erf(x) = 2ir 1 I 2 Jq e u2 du and a 2 is the linearly extrapolated r.m.s. 
mass fluctuation in a sphere of radius tq. The latter is given by top-hat 



8 



filtering of the power spectrum as 



<J 2 = ! 1 x 



°0 Y 1 



1 + zJ {l+zf 



sin krQ cos kr$ 
(kr ) 3 (kr ) 2 



2 



P{k)dk, (8) 



where ro is given by ^tttqP = M and where P(k) is the power spectrum. Al- 
though this approach has been criticized as too simplistic, numerical simula- 
tions (Efstathiou et al. 1988; Efstathiou & Rees 1988; Carlberg & Couchman 
1989) have shown that it describes the mass distribution of newly formed 
structures remarkably well. Making the standard assumption of a Gaussian 
density field, Blanchard et al. (1992) have argued that it is an accurate 
description at least in the low mass limit. Since we are interested only in 
extremely low masses such as 1O 6 M0, it appears to suffice for our purposes. 

We choose 5 C = 1.69, which is the linearly extrapolated overdensity at 
which a spherically symmetric perturbation has collapsed into a virialized 
object (Gott & Rees 1975). Letting f g denote the fraction of all baryons in 
galaxies of mass greater than M at z, this would imply that 



fa ~ 1 - erf 



1.69(1 + 2?) 
V2a (M) 



0) 



if no other forces than gravity were at work. However, it is commonly 
believed that galaxies correspond only to such objects that are able to cool 
(and fragment into stars) in a dynamical time or a Hubble time (Binney 
1977; Rees & Ostriker 1977; Silk 1977; White & Rees 1978). Hence the 
above value of f g should be interpreted only as an upper limit. 

A common assumption is that the first galaxies to form have a total 
(baryonic and dark) mass M c ~ 10 6 Mq, roughly the Jeans mass at recom- 
bination. Blanchard et al. (1992) examine the interplay between cooling 
and gravitational collapse in considerable detail, and conclude that the first 
galaxies to form have masses in the range 1O 7 M to 1O 8 M , their redshift 
distribution still being given by equation (||). To keep things simple we 
will assume that all early galaxies have the same mass M c and compare the 
results for M c = 2 x 10 6 Af o , 1O 8 M and lO n M . 

Let R(z; z*) denote the radius of a shell at z that was created at z* by a 
galaxy of baryonic mass Mj, = f2f,M c as in Section ||. Then the naive filling 
factor, the total bubble volume per unit volume of the universe, is 



4 u( ^3 Pb d fg(z*) , , , n3 [ r(z z*) 3 df g (z*) 

-ttR(z; z*Y—-f dz. t = 4>*(l+zY / — -f dz*, 

3 M b d{-z*) J z (1 + z^fl 1 d(-z*) 

(10) 
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where 

0* ps 1600/i- 1 M 5 " 2/5 (O b /0.06). 

Clearly nothing prohibits </> from exceeding unity. This means that nearby 
shells have encountered each other and that certain volumes are being counted 
more than once. If the locations of the bubbles are uncorrelated, then the 
fraction of the universe that will be in a bubble, the porosity, is given by 

P = 1 - e - *. 



If the early galaxies are clustered rather than Poisson-distributed, this 
value is an overestimate. For an extreme (and very unrealistic) example, if 
they would always come in clusters of size n and the clusters would be much 
smaller than the typical bubble size of 100 kpc, then it is easy to see that 
P ~ 1 — e~^/ n . For more realistic cases, simple analytic expressions for P 
are generally out of reach. Since we expect the clustering to be quite weak, 
we will use the Poisson assumption for simplicity. 

The uppermost panels of Figures |5| and|^ contain P(z) for various param- 
eter values, calculated numerically from equation (pi]) using the numerical 
solutions for r(z;z*). It is seen that the lower mass in Figure |(2x 10 6 Mq 
versus 1O 8 M ) gives higher filling factors, so that the expanding shells fill 
almost 100% of space by z = 5 for three of the four choices of / s (5). In ||, 
we see that almost 20% of the baryons must be in galaxies by z = 5 to 
achieve this. The greater efficiency of small galaxies is to be expected, since 
0* oc M~ 2 / 5 . Although some parameters still yield the desired P ~ 100% by 
z = 5 in Figure |6|, using present-day masses like M c = 10 11 Mq fails dismally 
(not plotted) for all choices of the other parameters. Roughly, the largest 
M c that works is 1O 8 M 

As can be seen, the dependence on fa (dashed versus solid lines) is rather 
weak. 

In order to calculate ctq{M c ) from the fluctuations observed on larger 
scales today, we need detailed knowledge of the power spectrum down to 
very small scales, something which is fraught with considerable uncertainty. 
For this reason, we have chosen to label the curves by the more physical 
parameter / s (5), the fraction of all baryons that have formed galaxies by 
z = 5. The four sets of curves correspond to fractions of 50%, 20%, 10% and 
1%. These percentages should be compared with observational estimates of 



metallicity, as will be discussed in Section 3.2 



The second column of Table [l] contains the values o"o(M c ) necessary to 
obtain various values of f g (5), calculated by inverting the error function 
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bcdmfi 


bcdmfi 


btiltedfi 


btiltedfi 


1% 


3.94 


4.8 


3.5 


2.5 


2.0 


10% 


6.18 


3.1 


2.2 


1.6 


1.3 


20% 


7.92 


2.4 


1.7 


1.2 


1.0 


50% 


15.06 


1.3 


0.9 


0.6 


0.5 



Table 1: Correspondence between various ways of normalizing the power 
spectrum 

in equation (0). The last four columns contain the bias factors necessary 
to yield this value of o"o(M c ) for two choices of power spectra (CDM and 
n=0.7 tilted CDM) and two choices of cutoff mass (M c = 2 x 1O 6 M and 
M c = 1O 8 M ). Thus b = 7/cio(M c ), where we define 7 to be the ratio 
between a at M c and a at 8/i _1 Mpc = Performing the integral (||) 

numerically with the CDM transfer function given by Bardeen et al. 1986 
(BBKS), h = 0.5, = 1, <C 1 and an n = 1 Harrison-Zel'dovich initial 
spectrum gives 7 rj 19.0 for M c = 2 x 1O 6 M and 7 rs 13.6 for M c = 10 s M . 
Using the CDM transfer function of Bond and Efstathiou (1984) instead 
gives 7 ~ 18.1 and 7 ~ 13.7, respectively. The BBKS transfer function 
is more applicable here since it includes the logarithmic dependence that 
becomes important for very low masses. The BBKS transfer function with 
a tilted (n=0.7) primordial spectrum yields the significantly lower values 
7 f« 9.71 and 7 Rj 7.93, respectively. 

Basically, Table [l| shows that any of our values of / 9 (5) become consistent 
with a feasible bias factor for some choice of power spectrum and cutoff 
mass. 

3.2 IGM enrichment 

These values of f g (5) should be compared with observational estimates of 
metallicity, since if the stars in these early mini-galaxies produce the same 
fractions of heavy elements as do conventional stars, then these percentages 
are directly linked to the fraction of currently observed metals that were 
made before z = 5. Some of the enriched shells may be observable as 
quasar absorption line systems, as intracluster gas, and, indirectly, as in the 
metallicities of old disk and halo stars. 

Observations of iron abundances in intracluster gas by HEAO-1, Exosat 
and Ginga (e.g. Mushotzky 1984; Hughes et al. 1988; Edge 1989; Hatsukade 
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1989) have shown that most clusters have abundances between 25% and 
50% of the solar value. Einstein observations have showed the presence of 
a large variety of other heavy elements in the intracluster gas (Lea et al. 
1982; Rothenflug et al. 1984). Most of this gas and some of these metals 
are believed to be "primordial", since the gas mass in clusters is typically 
several times greater than the observed stellar mass in the cluster galaxies 
(Blumenthal et al 1984; David et al. 1990; Arnaud et al. 1991). 

There are indications that the most of these heavy elements may have 
been produced as recently as around z = 2—3, and that the metallicity in the 
halo gas of some z ~ 3 galaxies inferred from QSO absorption line studies 
are as low as 0.1% of the solar value (Steidel 1990). However, this and other 
observations of extremely metal-poor objects (Pettini et al. 1990) does not 
necessarily rule out our scenario, since it is highly uncertain whether all 
the hydrogen in the swept-up IGM would get thoroughly mixed with the 
metal-rich supernova ejecta. 

3.3 IGM temperature 

Let T(z; z*) denote the temperature of the interior of a bubble at z that was 
created at z* as in Section ^. Then the volume-averaged temperature of the 
IGM is 

When <p becomes of order unity, the IGM swept up by the expanding shells 
is no longer cold, neutral and homogeneous, so the treatment in Section || 
breaks down. The resulting temperatures will be underestimated, since less 
thermal energy needs to be expended on heating and ionization. As can be 
seen in Figures || and ||, the transition from (f> <C 1, where the treatment 
in Section ^ is valid, to ^ > 1, where the IGM becomes fairly uniform, is 
quite rapid. Since Tjgm defined above is proportional to the thermal energy 
per unit volume, energy conservation leads us to assume that Tiqm remains 
fairly constant during this transition and therefore is a good estimate of the 
bulk IGM temperature immediately afterwards. From this time on, we will 
approximate the IGM outside the scattered dense and cold shell remnants 
by a uniform isothermal plasma. Applying equations (|2|), (^) and (||) to the 
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IGM yields the following equation for its thermal evolution: 

T 5 + e mj , (12) 



— T 

dz 5 



A(l + zf' 2 



l + z 



where the first term encompasses cooling from adiabatic expansion and the 
second term Compton cooling. T5 = Tjqm /10 5 K, A = l.5(to/t cornp ) « 
0.0042/i _1 and £i n j = t^Lmj/k x 10 5 iT, where Li n j is the power injected 
into the IGM per proton from all heat sources combined. The Compton 
cooling term is seen to increase with redshift, equaling the adiabatic term 
at x ~ 17h 0A - 1. 

In the most pessimistic case of no reheating whatsoever, i.e. for £i n j = 0, 
equation (|l~2^) has the solution 

To,(l + z) 2 e 0AA ^ 5/2 . 

A more optimistic assumption is that some fraction fi n j of the total energy 
released from stellar burning in newly formed galaxies continues to heat the 
IGM, i.e. 



/ 0.02 x 0.007m p c 2 \ df g ^ / 2.1xl0 4 \ / 1 

mj ~ hnj { k x 10^ ) d(-z) ~ I a i Jmje P { 2 



1.69(1+2) 



cr (M c 



which would also incorporate other modes of energy injection such as radia- 
tion. The middle panels in Figures || and || show the temperature resulting 
from fi n j = 0, 0.001 and 0.02 and different initial values. For P < 0.8, Tjgm 
has been calculated numerically from equation ([ll]) by using the numerical 
solutions for r(z; z*). Then the value at the redshift for which P = 0.8 has 
been used as initial data for equation (|l~2|). For comparison, two horizontal 
lines have been added showing what temperatures would be required to ob- 
tain neutral fractions of 10 -6 and 10~ 5 in equilibrium, using equation (|l3|). 
Since the plasma is in fact out of equilibrium, these highly ionized states can 
be maintained at much lower temperatures, as is seen in the bottom plots. 



3.4 IGM ionization and the Gunn-Peterson effect 

Assuming that any neutral hydrogen in the remains of the shells will have 
insufficient time to diffuse far into the hot ionized regions that used to be 
shell interiors, we can treat the latter as an isolated mixture of gas and 
plasma where the ionization fraction x evolves as 

X = n X [(1 - x)A ion (T) - X Kec(T)} , 
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and where the rates for collisional ionization and recombination are given 
by (Stebbins & Silk 1986) 



Aion = (crdv) « 7.2 al 



kT\ 1 ' 2 



m. 



-E /kT 



A, 



64tt 4 2 fkT\- 2 / 3 
:« a c I — I 



3\/37r u \EqJ 

where ao is the Bohr radius and a is the fine structure constant. Changing 
the independent variable to redshift, the $7 = 1 case leaves us with 



d 3 fmP-b 



dz^ 



2 (1 - 5) 3 

\on{T) = n c otoAion(7 1 ) 



(1 + z) 1 /2 x [(i _ x )A lon ( x ) - X A rec (T)], 



5.7 x 10 4 /iO f 



,1/2 1.58/Ts 



-2/3 



A rec (T) = n c0 toA rec (T) « [0.16M4] T 5 " 

For large enough z, the ionization fraction will adjust rapidly enough to 
remain in a quasistatic equilibrium and hence be given by \ = 0, i.e. 



X 



1 + 



ArecjT) 
Aion(T) 



1 + 2.8 x lo-eiT^e 1 - 58 /^ 



(13) 



The observed absence of a Gunn-Peterson trough in the spectra of high- 
redshift quasars strongly constrains the density of neutral hydrogen in the 
IGM. The most thorough study to date, involving eight quasars (Steidel & 
Sargent 1987), concluded that 

n Hl (z = 2.64) < (1.2 ± 3.1) x KT 8 /^ 1 

if O = 1. In our model this corresponds to (1 — x) < (1-2 ± 3.1) x 
10~ 8 /(/rA) ~ 2 x 10" 6 for U b = 0.06 and / m =0.1. Thus we are helped not 
only by the IGM being ionized, but also by it being diffuse. In a recent study 
of a single quasar, Webb et al. (1992) find the data consistent with either 
^/fj(-z = 4.1) = or 0,Hj(z = 4.1) = 1.5 x lO -8 /^ 1 , depending on model 
assumptions. We will use the latter value as an upper limit. Finally, recent 
Hubble Space Telescope spectroscopy of 3C 273 has been used to infer that 
$Ihi{z = 0.158) < 1.4 x 10 _7 /t^ . The constraints from these three studies 
are plotted in Figures || and |6] together with the ionization levels predicted 
by our scenario. 
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To achieve x = 10 4 ,10 5 and 10 6 in equilibrium would by equa- 
tion (ID require T > 5.5 x 10 4 K, T > 1.1 X 10 5 K and T > 3.6 x 10 5 iT, 
respectively. As can be seen from the numerical solutions in the bottom 
panels of Figures || and ^, the recombination rate is generally too slow for 
equilibrium to be established, and the IGM remains almost completely ion- 
ized even when T <C 15, OOO-RT and equilibrium would have yielded x ~ 0- 
In both 3a and 3b, a very moderate reheating (fi n j = 0.001, heavy lines) is 
seen to suffice to satisfy the three observational constraints. In the absence 
of any reheating whatsoever, the only models that satisfy the constraints 
are those with very low density = 0.01 or f m = 0.01). 

In summary, the only parameters that are strongly constrained by the 
Gunn-Peterson test are M c and ao(M c ). 

3.5 Other spectral constraints 

Let us estimate to what extent Compton cooling of the hot plasma will 
distort the cosmic microwave background radiation (CBR). Since for T e 3> 
T 7 the Comptonization y-parameter (Stebbins & Silk 1986) 

/■to kJ 1 

Vc= %rn e a t cdt 

Jt m e c z 

is linear in the plasma energy density ^| + |^ kT e n e at each fixed time 
(T e = Tjgm), all that counts is the spatially averaged thermal energy density 
at each redshift. Since the former is simply £f(z; z*) times the density of 
injected energy fgfsn^bpc 2 , the calculation reduces to mere energetics and 
we obtain 

poo l z \ rz, 

Vc = y* J \ / v 7 ! + z e t {z; z*)dz dz*, 
J z d{-z*) Jo 

where 

y* = \fsn^l^ « 9 x io- 7 ^. 

o TflpVj 

The current observational upper limit on y is 2.5 x 10~ 5 (Mather et al. 1994), 
so even if we take / ff (0) as high as 100% and make a gross overestimate of the 
integral by making all our galaxies as early as at = 30 and by replacing 
£t(z; z*) by its upper bound 60% for all z, z*, our y is below the observational 
limit by three orders of magnitude for = 1. 

Now let us estimate the optical depth of the IGM. It has long been 
known that reionization can cause a spatial smoothing of the microwave 
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background as CBR photons Thomson scatter off of free electrons. Since 
n e = XlGMnb, the optical depth for Thomson scattering, i.e. the number 
of mean free paths that a CBR photon has traversed when it reaches our 
detectors, is 

n= a t xiGMn b cdt = / y/l + zxiGMdz, 

where 

r t * = |-/ m O b ^^ « 0.07fi 6 /i. 

Let us evaluate the integral by making the approximation that xigm in- 
creases abruptly from to 1 at some redshift Zi on . Then even for Zi on as 
high as 30, Tt ~ 7.9Mlb/ m « 0.02 <C 1 for our fiducial parameter values 
h = 0.5, fife = 0.06 and f m = 0.1, so the probability that a given CBR pho- 
ton is never scattered at all is e -0 02 ~ 98%. Hence this scenario for late 
reionization will have only a very marginal smoothing effect on the CBR. If 
the shells are totally ionized as well, then the factor f m disappears from the 
expressions above which helps only slightly. Then Zi on = 15 would imply 
that 8% of the CBR would be spatially smoothed on scales of a few degrees. 



4 Discussion 

We have calculated the effects of supernova driven winds from early galaxies 
assuming a Press-Schechter model of galaxy formation and a CDM power 
spectrum. The calculations have shown that reionization by such winds can 
indeed explain the observed absence of a Gunn-Peterson effect if a number 
of conditions are satisfied: 

1. The masses of the first generation of galaxies must be very small, not 
greater than about 1O 8 M . 

2. There is enough power on these small scales to get at least 10% of the 
baryons in galaxies by z = 5. 

3. Except for the case where fife is as low as 0.01, there must be some 
reheating of the IGM after z = 5 to prevent the IGM from recombining 
beyond allowed levels. 

4. The commonly used thin-shell approximation for expanding bubbles 
must remain valid over cosmological timescales, with the mass fraction 
in the interior remaining much less than unity. 
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Whether 1) is satisfied or not depends crucially on the model for struc- 
ture formation. This scenario is consistent with a pure CDM model and 
some low-bias tilted CDM models, but not with top-down models like pure 
HDM. 

Observations of nearly solar abundances of heavy elements in intracluster 
gas have given some support for 2), which is roughly equivalent to requiring 
that at least 10% of the heavy elements in the universe be made before 
z = 5 (or whenever <j) 3> 1). As discussed in Section [T^, the observations of 
some extremely metal-poor objects in QSO absorption line studies do not 
necessarily rule out our scenario, since it is highly uncertain whether all 
the hydrogen in the swept-up IGM would get thoroughly mixed with the 
metal-rich supernova ejecta. The fact that large numbers of mini-galaxies 
are not seen today need not be a problem either. Possible explanations for 
this range from mechanisms for physically destroying them (Dekel & Silk 
1986, for instance) to the fact that the faint end of the luminosity function 
is still so poorly known that old dwarf galaxies in the field may be too faint 
to see by the present epoch (Binggeli et al. 1988). 

To violate 3), the actual reheating would have to be extremely small. A 
current IGM temperature between W^K and 10 5 K suffices, depending on 
other parameter values, since the low density IGM never has time to reach 
its equilibrium ionization. 

The thin-shell approximation 4) is obviously a weak point in the analysis, 
because of the simplistic treatment of the dense shell and its interface with 
the interior bubble. For instance, could the shell cool and fragment due to 
gravitational instability before it collides with other shells? An approximate 
analytic model for such instability has been provided by Ostriker & Cowie 
(1981). Their criterion is that instability sets in when S > 1, where 

„ _ 2Gp s heiiR 2 



Rv, 



and the sound speed v s = ^j5kT/3m p . In terms of our dimensionless vari- 
ables, this becomes 



0.067 vo.iy VI + z 



-1 / i , ,\3J 



0.011 x M^T- 172 fQb U 6 V ( l + Z y^ 



which indicates that with our standard parameter values, gravitational in- 
stability does not pose problems even with fairly low shell temperatures. 
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The reason that the shell density is not limited to four times the ambient 
IGM density is that the jump condition is not adiabatic, due mainly to 
effective Compton cooling at the high redshifts under consideration. 

After the critical z (typically between 20 and 5) at which the expanding 
shells have collided with neighbors and occupied most of space, the IGM is 
"frothy" on scales around 100 kpc, with dense cool shell remnants scattered 
in a hot thin and fairly uniform plasma. Since the dark matter distribu- 
tion is left almost unaffected by the expanding bubbles, formation of larger 
structures such as the galaxies we observe today should remain fairly unaf- 
fected as far as concerns gravitational instability. There is indirect influence, 
however: the ubiquitous metals created by the early mini-galaxies would en- 
hance the ability of the IGM to cool, which as mentioned in Section 3.1 is 
commonly believed to be crucial for galaxy formation. 

Blanchard et al. (1992) argue that if the IGM has a temperature higher 
than the virial temperature of a dark halo, pressure support will prevent it 
from falling into the potential well and thus stop it from forming a luminous 
galaxy. The virial temperature they estimate for an object of mass M formed 
at a redshift z is approximately 

K / M \ 2/3 
T m r « 5.7 x itf* {-^-) (l + z) 

for h = 0.5, so requiring T vir > Tjqm for say Tjqm = W 6 K at z = 5 would 
give a minimum galaxy mass of about 10 11 Mq. Such arguments indicate 
that the IGM reheating of our scenario might produce a "mass desert" be- 
tween the earliest mini-galaxies and the galaxies we see today: The first 
generation of galaxies, mini-galaxies with masses of perhaps 10 6 or 1O 8 M , 
would keep forming until their expanding bubbles had occupied most of 
space and altered the bulk properties of the IGM. After that, formation of 
galaxies much smaller than than those of today would be suppressed, since 
the IGM would be too hot. Eventually, as the IGM cools by adiabatic ex- 
pansion, a progressively larger fraction of the IGM can be accreted by dark 
matter potential wells. Indeed, even with the volume averaged IGM tem- 
perature remaining hot due to some form of reheating, cooling flows in deep 
potential wells, in particular galaxy clusters, would not be suppressed. Late 
formation of galaxies is therefore possible. 

Pressure balance between the shell and the interior during the expansion 
would give the ration T she u/T interior = Pb, interior / Pb, shell ~ 35/ m « 0.03 
for 5 = f m = 0.1, so the shell fragments would expected to contain non- 
negligible fractions of neutral hydrogen and thus absorb some Lyman-alpha. 
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A typical shell radius is about 100 kpc for M c = 10 M Q , a size comparable 
to that of the clouds of the Lyman-alpha forest. As to the number density of 
Lyman-alpha clouds, the observed velocity separations greatly exceed those 
we would expect if all shell fragments were to be identified with Lyman alpha 
clouds. Thus the majority of these fragments must have been destroyed by 
some other process. There are a number of ways in which this could occur, 
for instance through photoionization by UV flux from the parent galaxy or 
by collapse to form other dwarf galaxies. The resulting numbers resemble 
the abundance of minihalos in an alternative explanation of the Lyman- 
alpha forest (Rees 1986). Strong evolution, in the sense of an increasing 
cloud abundance with decreasing redshift, is expected to occur as cooling 
becomes effective. 

The authors would like to thank Alain Blanchard for discussions on the 
subject of the paper and David Schlegel, Douglas Scott and Charles Steidel 
for many useful comments. This research has been supported in part by a 
grant from the NSF. 
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Figure 1: Comoving radius of expanding shell. 
The comoving shell radius (1 + z)R is plotted for galaxies of total mass 
2 x 10 6 Mq, forming at integer redshifts from 1 to 29. Here $7 = 1, Qf, = 0.06, 
h = 0.5, and f m = 0.1. fd = 1 for the upper set of lines and fd = for the 
lower set. R has been truncated when T drops below 15,000 K, after which 
newly swept up IGM fails to become ionized. 
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Figure 2: This and the two following figures show the energy contents of an 
expanding bubble as a function of redshift, for different choices of fd and 
2*. CI = 1, Oj, = 0.06, h = 0.5 and f m = 0.1 for all three plots. Figures || 
and ^ illustrate the difference between fd = and fa = 1 (there is no shock 
cooling loss in the second case) . Figure || has fd = and illustrates that the 
Compton cooling loss is larger at higher redshift. 
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Figure 3: Energetics of expanding shell, example 2. 
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Figure 4: Energetics of expanding shell, example 3. 
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Figures [5] and ||: IGM evolution. 
Three different properties of the IGM (filling factor, temperature and ion- 
ization) are plotted as a function of redshift for different choices of M c , f g (5) 
and fd- Q = 1, Qb = 0.06, h = 0.5 and f m = 0.1 in the two previous figures, 
[| and |(| In all panels, the different families of curves correspond to different 
values of / 9 (5); 50%, 20%, 10% and 1% from left to right, with the rightmost 
cases being omitted where they fail dismally. In the porosity and temper- 
ature plots (the upper two panels of || and |||), dashed lines correspond to 
fd = 1 and solid ones to fd = 0, whereas only the pessimistic fd = case 
is plotted in the ionization plots (the lower third). In the temperature and 
ionization plots (the lower two panels), the three branches of each curve 
correspond to the three reheating scenarios: fi n j = 0, fi n j = 0.001 and 
finj = 0.02. 
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